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We study a current-biased 0-7T-0 Josephson junction made by high-T c superconductors, 
theoretically. When a length of the tt junction is large enough, this junction contains 
a vort ex- ant i vortex pair at both ends of the tt junction. Magnetic flux carried by 
the vortices is calculated using the sine-Gordon equation. The result shows that the 
magnetic flux of the vortices is suppressed to zero as the distance between the vortices 
is reduced. By applying an external current, the orientation of the vortices is reversed, 
and a voltage pulse is generated. The current needed for this transition and generated 
pulse energy are calculated. Macroscopic quantum tunneling (MQT) in this transition 
is also studied. The tunneling rate has been evaluated by an effective Hamiltonian with 
one degree of freedom. 

KEYWORDS: high-T c superconductor, n junction, sine-Gordon equation, half vortex, macroscopic quan- 
tum tunneling 



§1. Introduction 

Static and dynamical properties of long Josephson junctions (LJJs) has been studied by 
many authors for practical applications. Particularly, kink propagation in LJJs has been 
studied both theoretically and experimentally on the basis of the classical (damped) sine- 
Gordon equationS fl 

(ft xx - <t>tt - sin - a<p t + f3(j) xx t = /. (1) 

Here, a and (3 are damping parameters due to the quasiparticle tunneling loss and the 
surface loss, and / = I/Iq is an external current scaled by the critical current Iq, and the 
spatial variable x is normalized by the Josephson penetration length Aj. Recently, quantum 



effects of a kink have been discussed through interference effects,™ and macroscopic quantum 



In this paper, we consider a long Josephson junction including a so-called ir junction, which 
has negative critical current. In such junctions, the phase difference shows more complicated 
properties than traditional LJJs. Hence, we expect new phenomena characteristic of the 
sine-Gordon field with infinite degrees of freedom. 

The negative critical current has been discussed first in Josephson junctions, which have 
an insulator layer with magnetic impurities.!' !) More recently, it has been proposed that 
7r junctions can be made by unconventional superconductors with non-s-wave symmetry. 
Geshkenbein et a/.B have discussed it junctions formed in heavy electron systems, whereas 
Ti junctions in high-T c superconductors have been proposed to explain the positive param- 
agnetic Meissner effect.!'!'' In order to probe the symmetry of the superconducting gap, a 
number of experiments have been performed, including interference measurements in single 
crystal YE^CuaO^^-Pb SQUIDs,0'0 and direct imaging of magnetic flux in a tricrystal 
ring geometry.^ Ill) These measurements indicate that tt junctions can be realized at a grain 
boundary of high-T c superconductors, and that the (i-wave symmetry is realized in high-T c 
superconductors. 

Static behavior of LJJs with both and tt junctions (0-tt junctions) has been studied 
theoretically by several authors.SBIll) They have shown that a half vortex appears spon- 
taneously at a boundary between junction and tt junction in sufficiently long junctions. 
Further, Kuklov et al. have proposed that the half vortex can change its orientation by ap- 
plying an external current to the junction.© They have also conjectured that the half-vortex 
can be utilized in superconducting memory and logic devise. 

In this paper, we consider a 0-7T-0 long Josephson junction, which has a positive critical 
current at \x\ > a, and has a negative critical current at |x| < a. This junction may be 
realized in systems such as grain boundary Josephson junctions.© First, we study static 
properties of the junction in § |2|. After that, we consider MQT in this junction to study 
quantum effects in § |3|. Summary is given in § |j. 

§2. Hamiltonian and Static Properties 

2. 1 Hamiltonian 

The Hamiltonian of the 0-7T-0 Josephson junction with an external current is given by 



tunneling (MQT).|) 
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Fig. 1. Two types of solutions of (|4|) are drawn for a sufficiently large a. The minima of the potential for 
<j> is given as (2n + l)ir at \x\ < a, and as mr at |a;| > a. 



Here, x is scaled by Aj, and / = I/Iq is an external current, and the Hamiltonian is scaled 
by an energy scale of the junction, E . 0(x) is a step function defined by 

e(x) = ( - 1 (N<a) . (3) 

V ; \ 1 (\x\ > a) V ' 

Here, to simplify the situation, we assume that the absolute value of the critical current is 
uniform along the junction. Static behavior of is given by 5H/5(j)(x) = 0. From we 
obtain the sine-Gordon equation© 

4>xx = 6(x)sin0 + /. (4) 

This equation is also obtained from (|I|) by taking <ft t = 0. It should be noted that dissipation 
does not affect static properties. 



First, we consider the case with / = 0. In stable solutions by minimizing (0), the system 
prefers to have a uniform phase difference at 2irn for x 3> \a\, while = (2n + 1)tc is 
preferred for \x\ -C a, where n is an integer. As a result, we obtain two types of solutions, 
0i (x) and 02 (x) as schematically shown in Fig. 1. Here, we assumed without losing generality 
that both solutions satisfy 0(±oo) = 0. 

The physical meaning of these solutions is clear in case of a ^> 1. For example, we consider 
the situation described by 0i(x). The induced magnetic field in the junction is proportional 
to (p x . Hence, a vortex and an antivortex appear at x = —a and at x = a, respectively. The 
magnetic flux carried by these vortex is ±$o/2, where $o = h/2e is the unit flux. These 
induced vortices are called half vortices. The other solution 02 (x) has the same vortices, 
except that the orientations of the vortices are reversed from those in 0i(x). 
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Fig. 2. A<fi = (pi (0) — 4>i (00) is shown as a function of a, which is a half of the distance between two vortices 
scaled by Aj . The critical current f c necessary for vortices to change their orientation is also shown. Both 
varnishes at a — 7r/4 0.79 as a decreases. 

2.2 Magnetic flux 

The induced magnetic flux is proportional to A<f) = <fii(0) — <pi(oo). In the case of a ^> 1, 
we obtain A<p = ir, where the magnetic flux carried by the vortices is ±$ /2- The magnetic 
flux defined by A0, however, decreases as we put vortex and antivortex closer each other. 
We have calculated (f){x) from ([|) numerically, and obtained A0 as a function of a. We show 
the result in Fig. 2. For a < tt/A w 0.79, we obtained A0 = 0. This means that the vortices 
disappear when the distance between vortices is too small. As shown later, the critical value 
7r/4 can be obtained analytically. For a > 7r/4, the value of A0 increases quickly as the 
value of a becomes large, and for o>l, approaches ir. 

We expect that this behavior of the flux carried by the vortices may be applied to accurate 
measurements of Aj. When we change the length of tt junction or the Josephson penetration 
length Aj by external magnetic field, measurement of magnetic flux carried by vortices will 
tell us the information about the value of Aj. The measurement of Aj using direct imaging 
by scanning SQUIDs has already been reported by Kirtley et. al.® Compared with this, we 
expect that the measurement using the 0-7T-0 Josephson junction is more accurate, because 
the magnetic flux is sensitive to the ratio a = d/2Aj, where d is a length of the 7r junction. 

2.3 External current 

Next, we study static behavior of the solution in the presence of the external current /. We 
assume that the initial state is described by 4>i(x). As the external current / adiabatically 
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increases from zero, the form of <f>i(x) is modified. The modified solution is denoted as 
4>i(x, /). When the external current takes a critical value f c , the solution 0i(x, /) becomes 
unstable. Then, a transition from 4>i(x, f) to the other stable solution 2 (x, f) occurs. Here, 
4>2{x,f) is a solution modified adiabatically from the initial solution <p2{x) by the external 
current /. During the transition, a voltage pulse across the junction is generated. This 
transition is intuitively understood easily as follows: the vortices exchange their locations 
each other. After that, the system remains the state described by faix, f), as long as / > 0. 
When a negative external current is applied to the junction, a transition from <f>%(x, f) to 
0i (x, f) occurs at / = -f c . 

The critical current f c is calculated from (|J) numerically. The result is shown also in Fig. 2. 
For a > 7r/4, the critical current f c increases as a becomes large, and is saturated toward 
2/tt, which can be obtained analytically as shown later. Note that the critical value f c is 
smaller than 1. This means that the exchange of a vortex and an antivortex occurs before 
the whole junction is driven to a voltage state. 

2.4 Limiting cases 

Next, we consider two limiting cases: (i) a = tc/A + A (A 1), and (ii) a ^> 1. In these 
cases, we can perform analytical calculation. 

In the case (i), the phase difference <j)(x) satisfies |0(x)| <C 1 for all x. Then, the Hamilto- 
nian (||]) for / = up to <ft 2 is reduced to 

H[<j>{x)] = £\la^(*)i (~ + e(x)j <p{z), (5) 

where (p represents fluctuations around the trivial solution and is defined by <f>(x) = + (p(x). 
The eigenmodes of the fluctuation are obtained by solving the 'Schodinger equation' 

dx 2 

under the normalization condition 



+ Q(x)tp n = e n ip n , (6) 



dx ip n (x)ip m {x) = S nm . (7) 

Here, Q(x) can be regarded as a well-shaped potential with its width 2a. When a < 7r/4, 
the lowest energy Eq is positive, and the trivial solution = is stable. However, when 
a = 7r/4 + A(0<A^Cl), there exists a negative eigenmode 0o> which means the instability 
of the trivial solution <p(x) = 0. Hence, when we expand the static solution 

00 

(j)(x) = J2 C nfn(x), (8) 
n=0 
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the coefficients C n become for n > 1, and only Cq is nonzero. For a = 7r/4 + A(A^Cl), 
the ground state energy e is close to zero. By taking e = 0, the form of the unstable mode 
(fio(x) is obtained approximately from (|) and (JJ) as 



-cosx (\x\ < a), 

7T + 4 

' cosae^' 11 '^ (\x\ > a). 



(9) 



V 7T + 4 

To determine Cq, we derive an effective Hamiltonian for Cq by substituting <p{x) = Cq<Po(x) 
and (0) to the Hamiltonian (0). By assuming that Cq is small, a simple analysis gives 



H = -f c » + - V^fi /C ° (10) 

to the fourth order of Co. Here, £o = 8A/(7r + 4). From (|I0D, and / c are calculated 
analytically as 

= / c = -^-r— — rA 3//2 . (11) 

y V vr + 4 ' J 27(vr + 2) v ; 

In the case (ii), i.e. for a ^> 1, we can regard the junction as two independent 0-7r junctions. 

Then, a half vortex and a half antivortex appear at x — ±a. We notice only the half vortex 

at x — a. The behavior of a half vortex in 0-tt junctions in the presence of an external 

current has already been studied by Kuklov et.a/S They have studied a critical current / c , 

at which a half vortex changes its orientation with an integer flux being created. In order to 

study the stability of the static solution <fi\(x, f), we expand <f>(x) around <f>i(x, f) 

oo 

</>0r) = <Mx,/) + £C>„(x) (12) 

n=0 

with eigenmodes of small fluctuation defined by 

~ + u ' >iPn = En ^ u " = W iMx ' f)) - (13) 

Here, U(<f>) = Q(x)(l — cos0). At f — f c , the lowest energy Eq becomes zero as in the 
case (i). It can easily be checked that (p = Cd x (pi(x, f) always satisfies the equation fll~3|) 
with E — 0. Here, the constant C is determined by solving the normalization condition 
(0). This mode is called the translational mode, because this modulation translates the half 
vortex along the junction. We should note, however, that both <po(x) and d x (po(x) must be 
continuous at x — a. This condition and (|J) leads to 0i(a, / c ) = 0. Further, from the first 
integral of (|J), we obtain 

- cos 0(0) - f<P(0) = l -<p x (af - cos0(a) - /0(a), (14) 

+ cos0(oo) - /0(oo) = ^4> x {a) 2 + cos 0(a) - /0(a). (15) 
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Fig. 3. Sketches of the transition from <j>i(x) to <f>2{x) in two limiting cases: (a) a = ir/4 + A (A <C 1), and 
(b) for a» 1. The gray lines represent intermediate configurations during the transition. The transition 
occurs symmetrically in the case (a), while the integer flux is generated at either side in the case (b). 

From these equations and 4>\{a, f c ) = 0, the critical current for a ^> 1 is calculated as 
f c = 2/71.0* For the critical solution <j)(x, / c ), the normalization constant in <po is numerically 
calculated as C ~ 0.5468. 

As well as in the case (i), only the lowest energy mode (fo(x) becomes unstable for / > f c . 
The effective Hamiltonian in the case (ii) can be made by substituting (ft = <fii(x, f c )+C' ip (x) 
to @ as 

H = 7i5CC - a{CC ) 3 , (16) 

to the third order of C' . Here, 5 = f c — f denotes the difference between the critical current 
and the external current, and a is a constant calculated by 

1 f°° 

a = - dxQ(x)sin(j) 1 (x, f c )(d(f)i(x, f c )) 3 

6 J-oo 

= ^(0oo sin + cos 0oo - 1) w 0.1403, (17) 

where (poo = 0(oo) = — sin -1 f c . From the effective Hamiltonian (0), it is shown that there 
exists a metastable state for 5 > 0. Note that (|16D is valid only for a small C' Q . Hence, fllSD 
can be used only for 5 <^ 1 case, where the value of C at a metastable state is small. 

2.5 Crossover region 

It should be noted that the way to exchange the location of a vortex and an antivortex 
each other is different in the two limiting cases. In the case (i), the transition from <pi to 
02 occurs by keeping the symmetry <f)(—x) = 4>{x) as shown in Fig. 3(a), because the most 
unstable mode ipo(x) is symmetric. In the case (ii), the transition occurs symmetrically in 
the original model @ as well as in the case (i). However, in the presence of even small spatial 
inhomogeneities, either of two half vortices begins to move at a lower external current as 
shown in Fig. 3(b). Then, this half vortex changes its orientation first and creates an integer 
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Fig. 4. A sketch of the potential U"((f>i(x)) in eq. (|l8|). Here, E a is the ground state energy, and E\ is the 
first excited state energy. 

vortex in the region \x\ < a. This integer flux propagates along the tt junction toward the 
other half vortex, and combines to the other half vortex to change its orientation. 

To study the transition process for an intermediate value of a, we focus on the Schrodinger 
equation for fluctuations <p(x), 

- + u "<p» = E ^n, u" = f^OK*))- (is) 

The potential term U"(<f)(x)) has two minima at x = ±a as shown in Fig. 4. The ground 
state energy E and the first excited state energy Ei are also drawn in Fig. 4. As a increases, 
the wave function of the ground state is modified by keeping the symmetry (p (—x) = <fo(x) 
and its nodeless form. Hence, the most unstable mode <fo(x) is connected from the case (i) 
to the case (ii), smoothly. 

When the distance between two wells 2a increases, the energy splitting A = E\ — E is 
suppressed exponentially. For a>l, the lowest two eigenstates become almost degenerate, 
and we can constitute wave functions localized at each well as 

(pm = ( ( Po{x)±( Pl (x))/V2. (19) 

In the presence of spatial inhomogeneities, the potential energy is modified, and the energy 
difference between two wells appears. This effect can be studied by the effective Hamiltonian 
on the two-dimensional Hilbert space spanned by tpR and ipL' 

H cS = Aa x + ea z . (20) 

Here, e is the energy difference between wells due to inhomogeneities, and a x and o z are 
Pauli's matrices. When the energy splitting A is much larger than e, the wave function of 
the ground state is symmetric, and given by (fo approximately. Hence, the transition from 
<pi to 02 occurs symmetrically as shown in Fig. 3(a). On the other hand, A becomes smaller 



Fig. 5. The solid line represents AE — H[4>i(x, f c )] — H[(j>2(x, / c )] as a function of a. Here, Eq is an energy 
scale of the junction, and is defined below (|2l]). The curve for AE/E approaches 8a(gray line) for a ^> 1. 

than e when a is large. Then, the wave function of the ground state is localized in either well, 
and given by (pn (or ipi) approximately. As a result, the transition begins from either side 
as shown in Fig. 3(b). The qualitative change which is expected to occur in the intermediate 
region of a is not a transition but a crossover. 

2.6 Voltage pulse 

In order to estimate a voltage of the pulse, we have calculated numerically the energy 
difference AE defined as 

AE/E = H[Mx, fc)} ~ H[4> 2 (x, / c )], (21) 

where E = § I \j/27rL is an energy scale of the junction, and L is the junction length. The 
result is shown in Fig. 5. The energy difference increases monotonically for a > tt/4. When 
a is large, the curve for AE/Eq approaches 8a. This feature can be obtained analytically by 
neglecting the spatial distribution of vortices around the boundary between the junction 
and the it junction. The electric power of the pulse P is estimated as P = AE/T, where 
T is the time scale of the transition from X to 02 or vice versa. Accurate estimate of T is 
difficult, because we have to solve the sine-Gordon equation ([]]) dynamically. For a 1, 
however, it is inferred that T is determined by the time for an integer vortex to propagate 
from x = a to x = —a or vice versa, and is estimated as T ~ 2aAj/c in the original unit. 
Here, c = Aju; c is a characteristic velocity of the integer vortex, and u c ~ u v x max(l, a+/3/3) 
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is a characteristic frequency of the junction. Thus, the pulse voltage V is estimated as 

P $o^Aj 

7 Z — ' ^ ' 

where L is the junction length. Here, we used the critical current f c = I / 1$ = 2/tt. From 
Fig. 5, it is expected that the pulse voltage is suppressed as a decreases. Hence, the estimate 
(f22|) is expected to give an upper limit for V. 



§3. Macroscopic Quantum Tunneling 

In this section, we consider the transition due to macroscopic quantum tunneling (MQT) 
from the metastable state (j>\ to the stable state 02- in this paper, we neglect dissipation 
effects, taking a and (3 in (p]) as 0. Moreover, we study only limiting cases, which allow us 
to perform analytical calculation on the basis of the effective Hamiltonian (|T0|) and ( |16| ) . 

We first consider the case a = 7r/4 + A(A <C 1). We assume that 8 = f c — f <C 1, and that 
( |10"D can be approximated as a cubic potential. For a cubic potential, the tunneling rate V 
is calculated a 

7.2V 

T = Aexp(-B), B = - . (23) 

Tiujq 

Here, V is the energy barrier, and uj q is the frequency of small oscillations around the 
metastable state. The prefactor A is the order of u p , where u p is the plasma frequency of 
the junction. From (|T0D, we can estimate the exponent 

B « ^\-^5 5 /\ (24) 

where [3 2 = hu p /E . 

In a similar way, the exponent B is estimated for the case a ^> 1 from the effective 
Hamiltonian (fl6|) as 



B « (25) 



99 

J 2 

Note that the exponent B is proportional to 5 5//4 in both cases as seen in (p4|) and ( p5|) . 
Therefore, we expect that the ^-dependence in B does not change qualitatively in all the 
ranges of a. 

The value of ft 2 can be related to experimental parameters as 

1/2 

137 V W 2 e r 



?^{ 2 M\ • (^) 



where is the London length, and e T /d is the capacitance per area, and W is the junction 
width. The typical experimental value for j3 is very small (~ 10~ 3 ).& Therefore, MQT can be 
observed only for 5 1, because the tunneling rate V must be large enough to be observed 
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in the laboratory. For example, assuming a ^> 1, (5 = 10 3 and A = 10 10 [l/s], we obtain 
r~2x 10 2 [l/s] for 5 = 1(T 3 . 

The quantum effects were treated here on the basis of the effective Hamiltonians including 
only one degree of freedom for the field. Other degrees of freedom appears only in the 
form of plasmons in the semiclassical approximation, while the plasmons do not affect the 
tunneling rate at sufficiently low temperatures.B' The semiclassical approximation is justified 
for /3 < 1. Then, many-body effects characteristic of the sine-Gordon fields do not appear 
in this junction. 

In the above calculation, we have neglected dissipation effects due to quasiparticles. It 
is, however, expected that damping effect on the junction remains even at sufficiently low 
temperatures, and strongly affects the tunneling rate, when there exist gapless nodes for 
quasiparticle excitation in non-s-wave superconductors. Because dissipation on the junction 
generally suppresses the tunneling rate, our calculation gives an upper limit for Damping 
effects on MQT in this case with more accurate estimate of the tunneling rate remains for 
further studies. 

§4. Summary 

In summary, we studied static properties of magnetic fluxes and their macroscopic quantum 
tunneling in the 0-7T-0 Josephson junction, where two half vortices are formed if the tt junction 
region is long. We calculated the magnetic flux of spontaneously induced vortices, and the 
critical current needed to make a transition between two degenerate vortex configurations. 
We also studied quantum tunneling rate for this transition. This MQT may be observed in 
high-T c superconductors under an appropriate condition. 
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